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Abstract 

This paper provides a small noise approximation for local random center manifolds of a class 
of stochastic dynamical systems in Euclidean space. An example is presented to illustrate the 
method. 
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1 Introduction 

Center manifolds, together with stable and unstable manifolds, provide geometric structures that 
help understand nonlinear and stochastic dynamics. For stochastic dynamical systems, it is difficult 
to visualize or depict center manifolds. In the paper, we derive a small noise approximation for 
random center manifolds. 

'This work was done while Jian Ren was visiting the Institute for Pure and Apphed Mathematics (IPAM), Los 
Angeles, USA. It was partly supported by the NSF Grant 1025422, the NSFC grants 10971225 and 11028102, and 
the Fundamental Research Funds for the Central Universities (HUST No.2010ZD037). 
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We consider a system of nonlinear random differential equations (RDE) on a probability space 

(ill 

— =A{etUj)u + F{etu;,u), (1) 

where u € M", and 0^ : — ?> is an ergodic flow that preserves the probability measure P. Assume 
that is a fixed point, i.e., F(9tUJ,0) = for all t and all uj. 

This RDE system is usually coming from a stochastic differential equation (SDE) system, by a 
stationary coordinate transform [.5, 7]. 

Boxler [2, 3] proved a random center manifold theorem, and also considered reduction to the 
center manifold. Power series, in terms of variables on the center manifold, approximations are 
also available in [.'], S]. Sun [7] provided a small noise approximation for random stable or unstable 
manifolds. 

In this present paper, we derive a small noise approximation of random center manifolds, for 
the stochastic system (1) where the linear operator A{6tU}) has exponential trichotomy property. 



2 Exponential trichotomy 

Let (fi, J^, P) be a probability space with an ergodic flow 9t- A measurable map (j){t, ijJ,u) -.MxQx 
M"" R" is called a random dynamical system (RDS) if it satisfies the following conditions 

(i) 4>{0,uj,u) = u; 

(ii) (f>{t + s, u, u) = (p{t, 9s00, u)(j){s, uj, u), for all t, s € M and all a; € 0. 
The latter is the so-called cocycle property. 

For a linear RDE system in M" of the form 

f = A{9tu)u, (2) 

where matrix A{9tuj) satisfies the integrability condition that A € L^(0, J-", P). The integrability of 
A{9tuj) implies that the RDE system (2) generates a linear cocycle <I>(t, uj), which is the fundamental 

matrix e^^o ^(^'r'^)dr_ 'pj^gj^ ^j^g multiplicative ergodic theorem (MET), in [1, Chapter 3], yields that 
we have Lyapunov exponents Ai > A2 > • • • > A,-, with r < n, together with r Osledets subspaces 
Ei{uj), so that M" decomposes into the direct sum: 

= Ei{uj)®E2{uj)® ...Er{uj). (3) 

All the subspaces Ei{uj) are measurable and are random invariant subspaces of ^{t,uj), i.e. for 
i = 1,2,--- ,r, 

^{t,uj)Ei{uj) = Ei{9tuj), forallteM, F - a.s. (4) 
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We define stable, center and unstable subspaces 



E'iu) ^ (Bx,<oEi{u;), (5) 
^^=(0;) ^ ®x,=oE,{oo), (6) 
E^ioj) ^ ®x^^oEiioj). (7) 



Therefore, 

and ^{t,uj) can be decomposed as 

<^{t,uj) = $,(t,a;) e^>c(t,w) e$„(t,w). 



Denote 



then 



Xs = max Xi < 0, = min Aj > 0, 

A,<0 Ai>0 



lim - log \^{t,uj)v\ = X, 

f— >-±oo t 



with 

A = A„ for V e E'{io)\{0}, 

X = Xu, for V G ^"(tj) \{0}, 

A = 0, for V G E^{io) \ {0}. 

Taking a positive number 7 < ^ min{|As|, A^}, the MET leads to the following estimation [2, Lemma 
4.2]: 

\^s{t,uj)v\ < e(^''+^)*i;, v€E', t > 0, 
\<^c{t,uj)v\ < e'^l*!?;, v e E"", t G M, 

< e(^"-^)*u, uG-E", t<0. (8) 



Similar to the definition of exponential dichotomy in [•'">], the cocycle ^{t, u) with the conditions 
(3), (4) and (8) hold is said to satisfy the exponential trichotomy. 

For the nonlinear RDE system 

(111 

— =A{etu)u + F{etuj,u), (9) 
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we assume that the nonhnear term F is Lipschitz continuous on M", that is, 

\F{ui) - F{u2)\ < hp \ui - U2I, 

with the Lipschitz constant hp > 0. In the next section, we consider the random center manifold 
for this system. 

If the nonhnear term F is only locally Lipschitz, then we let F^^\u) = where 
Xr{u) is a cut-off function. Thus F^^^ is globally Lipschitz with Lipschitz constant Rhp; see [4, 
Lemma 4.1]. In this case we would obtain a local random center manifold. 

By using the fundamental function 
the solution of the random differential system (9) can be interpreted as 

u{t,UJ,Uo) = ^{t,Uj)uo + I ^{t - T,Uj)F{6rUJ,u{T))dT . 

Jo 

This equation has a unique measurable solution due to the Lipschitz continuity of F, and the solu- 
tion mapping {t,u},uo) 1-^ u{t,uj,uo) generates a random dynamical system (j){t,uj,uo). Moreover, 
the Jacobian mapping DuQu{t,LJ,UQ) = ^{t,uj) satisfies the exponential trichotomy. 



3 Random center manifold 

In this section, we consider random center manifolds for a SDE system with small multiplicative 
noise. We will recall a random center manifold theorem [2] and then approximate the random 
center manifold under small noise, in the case of non-positive Lyapunov exponents. This is done 
in deterministic case in [4] and [6, Chapter 3]. 

We introduce the definition of random center manifolds. 

Definition 1 (Random Center Manifold). A random set M{oj) is called a random center manifold 
for a random dynamical system (p(t,uj,x), if it satisfies the following conditions 

(i) It is an invariant set, i.e. (j){t,Lj, M(uj)) C M{6tuj) for all t. 

(a) It can be represented as a graph of a mapping from the center suhspace to its complement, 
i.e. there is a mapping h'^{uj, •) : — )■ E"^ © i?*, such that M{uj) = | (f , h'^{uj, -y)) : v € -E'^}, where 
h^{uj,0) = 0, h^{-,v) is measurable for every v ^ E^ and the tangency condition Dh^{uj,0) = 
holds. The center manifold M{uj) is often denoted as ]VF{uj). 

It is called a lipschitz center manifold if the mapping h'^{u;,-) : E^ E" is lipschitz and 

the tangency condition is absent. We call ]VP[uj) a local center manifold if it is a graph of the 
mapping Xr{v) h^{uj,-), where Xr{v) is a cut-off function. 
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For the nonlinear SDE system defined in the previous section, Boxler [2, •■]] has shown that there 
exists a (local) random center manifold, as the graph of a mapping that satisfies the following 
Liapunov-Perron equation 

h''{uj,x'') = (^j^ $;i(T,a;)F"(0^a;, ((/)^(r,^J,x^/^^(w,x'=)),/^^(0^a;,0^(r,w,x^/^'=(^^,2;^)))))dr, 

^ «>~^(r, tj)F" (Oruj, ((^"(r, w, h^{oj, x")), h^iOrUJ, 4>^{t, uj, x^, h^{uj, x^)))) ) dr^ , 

(10) 

for € and uj ^ il. where F*, are respectively the projection of F to the stable and unstable 
subspaces, (p^ is the projection of (j) to the center subspace. 



3.1 Center manifold of a system with small noise 

We consider a SDE system in Stratonovich form 

± = A''x + f''{x,y) + {ex^ oWiy, x G i?", 

y = A^y + nx,y) + {ey^oW,y, yeR^, ^^^^ 

where indicates the transpose of a vector or a matrix. In this system, A'^ and A^ are respectively 
n X n and m x m matrixes. The nonlinear functions f'^{x,y) : M" x — >■ and f^{x,y) : 
M"- X R"* are C"^-smooth and Lipschitz with Lipschitz constant Lf, i.e. f^{x,y) and f^{x,y) 

satisfy 

\nxi,yi) - nx2,y2)\ < Lf (|xi - X2I + \yi - ?/2|), 

and 

\f(.xi,yi) - f{x2,y2)\ < Lf (|xi - X2I + \yi - y2\), 

where |x| is the Euclidean norm of x G R" or x G R™. Usually, f^,/^ are locally Lipschitz 
continuous and in that case we would get a local random center manifold. The noise intensity e 
is a small positive parameter. Moreover, {Wf : t G R, i = 1,2} are n x n and mxm matrixes, 
respectively, with two-sided scalar Wiener process Wt as principal diagonal elements and all other 
elements being zero. 

For matrixes A^ and A*, we make the following assumption. 

(H): There are positive constants /3, 7 and satisfying /3 > 7 > and -ftT > 0, such that for 
every x G and y G i?™, the following exponential estimates hold 

|e^"*2;|Rn < Ke^\*\\x\Rn, t G R; |e^'*y|ij™ < i^e"^*|y|Km, t > 0. 

To facilitate random dynamical systems approach, we convert (11) into a system of random 
differential equations (RDE). To this end, consider linear stochastic differential equations 

dz, + Zidt = dWi, 1 = 1,2. (12) 
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Each of these equations has a unique stationary solution 

Zi{u:)= f em\, i = l,2. 

Moreover, 

Zi{etuj)= [ e^~'dW^, i = l,2. 

J — oo 

Thus, the hnear differential equations, 

dZi + Zidt = edWi, i = 1,2. (13) 

have unique stationary solutions Zi{uj) = £Zi{uj) and ZiiOfUj) = ezi{9tOj) for i = 1,2. 

Note that Zi{uj), Zi{6tuj), Zi{uj) and Zi{9tuj) for i = 1,2 are all principals diagonals matrixes, 
with the principal diagonals values e"^ dWV, j'i^ e"^"* dWr, £ dWV and e e"^"* dWV 

respectively. 

We introduce a random transformation 



Under this transformation, the SDE system (11) is converted into a RDE system 



\X = A^X + eziX + F%X,Y), XeR^, 
\y = A'Y + £Z2Y + F''{X,Y), YeR"". 

with 
and 

being Lipschitz functions with the Lipschitz constant Lf. 

For every 77 > 0, we define a product Banach Space := C,J x C^, where 

= {v? : M ^ M" : is continuous and sup | e~''"'l*le~=i'o ^i(^'-'^)dr-(^(t)|]g„ < 00}, 

tm 

C2 = {v? : R ^ : </. is continuous and sup | e-'^'I'le^^/o ^i(^'-'^)d'^yj(t)|jj„ < 00}, 
with norm 

\^{t)\ci = sup I e-'^l*! /o -iC^--) d-^(t) , 



(14) 
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and 

\cl){t)\c2 = sup|e-''l*le-^i'o^i(^-'")'iXt)|M.n. 

respectively. Furthermore, the norm in is 

\{X,Y)\c, = \X\ci + \Y\c2. 

The integral form of the differential equation (14) can be written as 
\Y{t) = f_^e'^°^^-^^+^-fr'^i'^rUj)drj^s^x,Y)dT. 



(15) 



Due to the Lipschitz condition of F'^ and in (14) and the measurability of Zi{uj), i = 1,2. the 
RDE system (14) has a unique solution. Therefore, the solution mapping generates a RDS. By 
the assumption (H) on A^, and the propersties of Zi{u), i = 1,2, as in [5], the exponential 
trichotomy holds (with E'^ empty). We thus have following result [2, 3]. 

Theorem 2. Assume that H holds and that j < rj < 13 satisfy the gap condition + < 1. 

Then there exists a random center manifold ^{u}) = (^, //^(cj, ^)) for (14), where H^{lo,-) : M" 
M™ is a Lipschitz function. 

3.2 Small noise approximation of random center manifolds 

We now approximate the above random center manifold under small noise intensity e. Denote the 
corresponding deterministic center manifold (when e = 0) and the random center manifold as 

H\i) = {{i,H\0)}, (16) 

and 

n^iv^i) = {{i,H%oj,i)) = {i,H\i) + eH\u;,i) + Oie'))}, (17) 

where H'^{-) and n(uj, •) are locally defined on M". Denoting x{d) = ^ € M", as in [5], the center 
manifold T-L^{uj,5,) can be converted back to the center manifold for the original stochastic system 

(11), 

= {(e"i('^)|,e"2M^^(a;,f))} = {(e,e"^2M^e(^^g-£^i(a;)^)^| 

= { (C, + eH\u,, + ez2{oj)H\0 - eHf{i)zi{u:)i + 0{e^))]. 

(18) 
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To approximate the center manifold, we expand 

X{t) = Xo{t) + eXi{t) + e^X2{t) + ■■■ , with X{0) = i = Xo{0), i G M", 

and 



with 



Noting that 



Y{t) = Yo{t) + sYiit) + e^Y2{t) + ■■■ , (19) 
y(0)= H'{u;,i)= H^{i)+ eH\u,0+ 0{e^). 



^-ezi = 1 _ + _^ + . . . ^ for i = 1, 2, 
we have, 

F%X{t),Y{t)) = e-^^i(^*'^)r(e"i(^*^)x(t),e^^2{^?*^)y(i)) 

= r{Xo{t),Yo{t)) + e[f^{Xo{t),Yo{t))[Xi{t) + zi{etUj)Xo{t)] 

+4^(Xo(t),FoW) [Yi{t) + Z2{etOj)Yo{t)] - ziietu;)r{Xoit),Yoit))} + 0{e^), 



and 

F'{X{t),Y{t)) = e-^'^^^"^^f{e'''^^'^^X{t),e^'^^^*'^^Y{t)) 

= r{Xo{t),Yo{t)) + e[f^{Xo{t),Yo{t))[Xi{t) + zi{etco)Xo{t)] 

+f;{Xoit),Yo{t)) [Y^it) + Z2ietu;)Yo{t)] - Z2(Mr (^o(t), ^oW) | + 0(8^). 



Inserting these expansions into (14) and matching the terms with the same powers of e, we get at 
the 0-th order 

(Xo{t)) ^ (A- Q\ (Xo{t)\ ^ (r {Xo{t),Yo{t))] 



\Yoit)J VO \Mt)J ^ \r{Xo{t),Yoit))J^ 

which can be expressed as 



(20) 
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Furthermore, at the first order in e, 



A" 
A' 



+ 



f^{Xo{t),Yo{t)) f^{Xo{t),Yo{t))\ 
f^{Xoit),Yo{t)) f^{Xoit),Yoit))J 



Xi 
Yi 



+ 



r^zi{etOj)Xo{t) + zii9tOj)XQit) + f;;z2{etu;)Yoit) - zii9tu;)r\ 
r^zi{etu)Xo{t) + Z2{etu:)Yo{t) + Z2{etu;)Yo{t) - Z2{etu;)rJ ' 



(21) 



which can be rewritten as 



Xi{t) 
Yi{t) 



exp 





A' 



Wf^{Xois),Yois)) f^{Xo{s),Yo{s)) 
\f^{Xois),Yo{s)) f^{Xois),Yois)) 



ds 



\ym) 



+ / exp 





A' 



(t-r) 



\n fy) J V/x^l^0 + ^2>0 + 4^^2>0-^2/V 



(22) 



Hence, 



J — oo 

/ e-^'"(l + e / Z2{eruj)dr + 

J — oo ^ 



e^{f:;'z2{9ruj)drf 
2! 



+ ••• {r(Xo(r),yo(T)) 



/,^(Xo(r),yo(r))(Xi(T) + zi(0,o.)Xo(t)) +4^(Xo(r),yo(r))(yi(r) + Z2(^.^)l"o(r)) 
Z2{erU;)r{Xoir),Yoir))] } dr + ©(e^) 



e-^'-r{XoiT),Yo{T))dT + e 



Z2{9rCo)drf{Xo{T),Yo{T)) 



+f^{Xo{T),Yo{T)) (Xi(r) + zi{erU;)Xo{T)) + 4^(Xo(r),yo(r)) (yi(r) + Z2(^ra;)lo(r)) 
-^2(^rO;)r(Xo(r),yo(r))] dr} + 0(e2). 

Therefore, matching with (19), we get 



H\0= f e-^V(^o(r),lo(r))dr, 



(23) 
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and 

H\uj,0 = f e-^^^{ /'%2(Mdrr(Xo(T),yo(r)) + /^(Xo(r),yo(r))[Xi(r) + zi(0,a;)Xo(r)] 

J —OO J T 

+4^(Xo(r), Yo{t)) [Yiir) + Z2(e.u;)yo(r)] - ^2(^.0;)/^ (Xo(r), Yoir)) } dr. (24) 

Moreover, by (18), the center manifold for the original stochastic system (11) is 
n^ico) = { (e, H%oj, 0) } = { (e, H\C) + sHHoj, + Oie')) } , 

with 

H'^iO = H'^iO, (25) 

and 

H\u;,0 = H\u:,0+Z2{uj)H\0-Hl{0^i{^)i- (26) 

We summarize the above approximation result in the following theorem. 
Theorem 3 (Approximation of a random center manifold). Assume that (H) holds and that r] sat- 
isfies J < rj < 13 and + < 1. Then there exists a center manifold 7i^{oj) = { (^, H'^{lo, ,^)) } 

for the stochastic system (11), where H" {uj , = H'^{i)+£H^{uJ,C)+0{e'^), withH'^{C) andH^{uj,^) 
in (25) and (26), respectively. 

We now look at an example of the approximation of random center manifolds. 
Example 1. Consider a SDE system 



X = a'^x + ex o Wt, x G M, 
ij = -y - x"^ + ey oWt, ye 



2 , _ ... . . (27) 



Here A" = a", A' = -1, f%x,y) = 0, f%x,y) = -x^, z{uj) = dVF^ and ziOtio) = 

e"^"* dWV = e~^z{uj) + e~* Jq e'^dWV- We restrict this system on a bounded disk contain- 
ing the origin. The new system satisfies the assumption (H). In this way, we obtain a local center 
manifold and we now consider its small noise approximation. 

We transform this SDE system into a RDE system, 
(xit) = a^X + ez{etu)X, x G M, 

[y(t) = -y + e2(0tt^)y - e-^^(^*^)(e^^(^*"')X)2, y G M. 
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Denote a = a'^ + £z{9tuj), b = —1 + ez{9tuj), and A = || ) . By the property 



^ lim J Jq z{6rOj) dr = in [5], we have the Lyapunov exponents 



\r = hm T log le-^o ^('^) 'i'^til = a'^ = 0, when u 



Ui 



t— S-±CXD 



As = — 1 < 0, when u 



(29) 



By (20) and (23), we get 



Xo(t) = H\i) = / e^-X^) dr = 



and 



Yo{t) = e-'H\l) - ~e{l - e-*) = -'f . 



From (22) and (24), we obtain 



Jo 



and 



H\uj,i) = f /%(0,,L^)dr(-X2(T)) -2Xo(T)(Xi(r) + z(0,L^)Xo(T)) +z(0,^^)X2(T)]dr 



OO J T 



e [ / ziQr^) dr + 2( / ziQr^) dr + ziQ^bj)) - z{dr^)\ dr, 
-OO Jt Jo 



thus 



H\u;,0 = T^Af f z{eri^)drdT- f e-z(0,a;) dr] = 0. 

J- + 2^ J -OO Jt J -OO 



(30) 



By (25) and (26) 
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and 



H\u,i)+z{u:)H\i)-Hl{i)z{u)i 
- z{Lo)i'^ + 2z{uj)f 



— oo 



Therefore, we have an approximation for center manifold of (27) for e sufficiently small 



with 




■0 



— oo 



The center manifold of the corresponding deterministic system 



X = 0, 

y = -y- x^- 



(31) 



is the graph of the function /i'^(^) = A few realizations of the approximate center manifold, 

together with the deterministic center manifold, are shown in Figure 1. 
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